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STRANGE STARS : 
WHICH IS THE GROUND STATE OF QCD 
AT FINITE BARYON NUMBER ? 

Dany PAGE ^ ^ 
Astronomy Department, Columbia University 
New York, NY 10027, USA 

Abstract 

Witten's conjecture about strange quark matter ('Strange Matter') being 
the ground state of QCD at finite baryon number is presented and stars made 
of strange matter ('Strange Stars') are compared to neutron stars. The only 
observable way in which a strange star differs from a neutron star is in its 
early thermal history and a detailed study of strange star cooling is reported 
and compared to neutron star cooling. One concludes that future detection 
of thermal radiation from the compact object produced in the core collapse 
of SN 1987A could present the first evidence for strange matter. 

1 Introduction 

In a seminal 1984 paper E. Witten^ proposed that the ground state of baryon 
matter (or QCD) at finite baryon number might be quark matter and not nuclear 
matter as always supposed. If the strange quark mass is neglected it is a simple 
consequence of the Pauli principle that three fiavor quark matter has a lower energy 
than two fiavor quark matter. The conclusion usually still holds when is taken 
into account. (Eour fiavour is forbidden due to the much larger mass of the charm 
quark). Quark matter must then be three fiavoured, and is called Strange Matter. 
This simple fact was of course well known long before Witten's paper. What was 
new is that strange matter may have a lower energy per baryon than nucleon matter. 
In light of this on can phrase Witten's conjecture the following way : 
Two flavor quark matter, at zero pressure, has a higher energy per baryon 

than the ground state of nucleonic matter (Fe^^ with E/A = 930 MeV) 
but three flavor quark matter ('strange matter') has a lower energy per 
baryon. 

Eor small A Witten's hypothesis is certainly wrong, e.g. for A = 1 we obtain 

for 'strange matter' the A which has a mass of 1115 Mev, which is higher than 
the nucleon mass by quite a large amount. Eor slightly larger A's one can expect 
nucleonic matter to remain the ground state until some critical value. Eor heavier 
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nuclei it is easy to understand why they do not decay to the 'true' ground state 
: one must transform them to strange matter and not to two flavor quark matter, 
and thus must transform simultaneously A quarks u or d into s, an A^'^-order weak 
process. Nuclei are therefore only meta-stable but their lifetime far exceeding the 
age of the universe. 

Witten originally proposed strange matter as a candidate for dark matter to 
be formed in the early universe, but it was then shown that strange matter would 
not have survived very long and probably would not have formed at all at this 
stage. Presently, the most probable place where strange matter is expected to 
exist is in collapsed stars ('neutron stars'). Relativistic heavy ions collisions are 
also expected to produce quark matter, but not cold (degenerate) quark matter as 
we are discussing here. We refer to the review by Alcock & Olinto^ for a general 
discussion of the subject. 

We show here how a star made of strange matter, a 'strange star', can be distin- 
guished from an ordinary neutron star, and will conclude that the 1987A supernova 
could provide us with a unique oportunity to verify Witten's conjecture. 

2 A Model for Strange Matter 

In order to study quantitatively the properties of strange matter Farhi & Jaffe"^ 
have proposed a model based on the three following assumptions : 

1. The system is well approximated by a degenerate Fermi gas, separated from 
the vacuum by a phase boundary. The quark phase vacuum thus carry an 
energy density B (the 'bag constant') of the order of [lOOMeVy — (200Meu)'* 
(A typical bag model value is B = {14:5MeVy = blMeV/ fm^ ). 

2. Quarks are characterized by their current masses : 

mu = md = (1) 

m, ~ 100 - SOOMeV (2) 
(The prefered value is = ITbMeV). 

3. Interactions are modelled using renormalization-group-improved flrst order 
QCD, with a QCD coupling constant etc = qII^t^- 

The model thus depends on three parameters : B, nis and etc- Unfortunately neither 
of them is yet determined in the particular conditions considered here. However it 
allows us to test Witten's conjecture within a simple framework. 

The calculation starts with the grand thermodynamic potential (Landau's po- 
tential) fl(/ii)(= — Py), which has been calculated by several authors as a function 
of the various chemical potentials jii. The number densities are given as usual by 
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rii = —dQ/d/ii, i.e. rii = ni[iii). It is found that /i„ > Pf{u) and /id > Pf{d) 
{pf{i) is the Fermi momentum of the i particle), and that at not too high density 
/is < {ml , reflecting the fact that one gluon exchange is repulsive for 

relativistic particles and attractive for non-relativistic ones. Thus interactions favor 
the appearance of strange quarks. To find the chemical composition of the system 
at a given density ub one must solve the following system of four equations for the 
four jii [i = M, J, s, e) : 

flu + tie = f'd (3) 
tl's = f'd (4) 

2 11 

-riu- -rid- -ns-ne = {) (5) 

\^i^u + nd + n,s) = Ub (6) 

The first two equations impose beta equilibrium, the third charge conservation and 
the fourth one baryon number conservation. One still must impose the condition of 
zero pressure : 

P = -n„ -Qd-^s-^e- B = (7) 

which fixes ub- So, for each set of parameters (5, and etc) one obtains a chemical 
composition and an energy density e = Y^[Q,i-\-iiini)-\-B^ and so an energy per baryon 
E/A = e/uB which is compared to the value for nuclear matter. The result of the 
study is that strange matter does have a lower energy per baryon than nucleonic 
matter for quite large ranges of parameter values (see Fig.l in Ref."^). Of course the 
actual values of the parameters can perfectly well lie outside these specific ranges 
and the model used is questionable, but it makes Witten's idea plausible. The proof 
or refutation of Witten's conjecture can only come from experiment or observation. 

3 Strange Stars 

Since the most likely place to find strange matter is in collapsed objects ('neutron 
stars'), the study of strange stars has received quite a lot of attention, in the hope of 
finding some characteristic feature which would sharply differentiate a strange star 
from a neutron star. The basic properties of strange stars were first described by 
Alcock, Farhi & Olinto,'* and we sumarize here the most relevant ones. The structure 
of the star, assumed to be spherically symmetric, is obtained by integrating the 
Oppenheimer-Volkoff equation with the equation of state (F.O.S.), i.e. the relation 
between e, ub and P, given by the solution of equations 3 to 6 and equation 7 with 
P ^ 0. Unfortunately their global properties, such as mass and radius, turn out to 
be very similar to those of neutron stars and cannot be used to distinguish them. 
The (erroneous) announcement of the observation of a very fast rotating object in 
the remnant of SN1987A with a period of half a millisecond stimulated work on 
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rotational properties of neutron and strange stars, raising the possibility of putting 
strong constraints on the E.O.S. of hyperdense matter. However the observation 
was later rejected, ruining a second conceivable discriminant. 

We now review the basic structure of neutron stars and strange stars and show 
that they differ significantly in the crust. 

3.1 Neutron Star Structure : Crust 

A typical neutron star has a mass of f.4 Mq{Mq= 2-f0"^"^(/m is the sun mass) and 
a radius of the order of f km. Its average density is between 10^"* and 10^^ gm / cm^ ^ 
and as a zeroth approximation it can be seen as a giant nucleus of ~ 10^'' nucleons 
and electrons (and muons) inside it to keep charge neutrality. At closer look it has 
a crust (mostly solid) made of nuclei on top of the nuclear matter. The uppermost 
layers are made of material accreted from the interstellar medium, but below this the 
matter has been processed to nuclear equilibrium, i.e. iron nuclei Fe^^. At densities 
above p ~ gm / cm^ the iron is fully (pressure) ionized and above p ~ 10^ gm/cm^ 
the electrons become relativistic. This latter fact implies that the (degenerate) 
electron gas contains a substantial amount of energy and when the density increases 
still further the electrons have to be absorbed into the nuclei, changing the chemical 
composition toward more neutron rich nuclei. At a density just below pdrip = 
4.3 • 10^^ gm/cm^ the equilibrium nucleus is Kr^^^ [N = 82, Z = 36), but above 
this point the neutrons produced cannot be all bound into the nuclei and start 
dripping out : we have now a material of neutron rich-nuclei immersed into a sea 
of neutrons and electrons. At higher densities the nuclei are deformed, becoming 
ellipsoidal, more and more elongated, form foils and then start merging together, 
giving a Swiss cheese-like structure, which dissolves into the nuclear soup around 
pnud = 2.8 • 10^'^gm/cm^. We thus have the outer crust at densities below pdrip^ then 
the inner crust in the neutron drip regime and below them the core at densities 
larger than pnud- 

3.2 Strange Star Structure : Crust 

Since strange matter is assumed to be the ground state of baronic matter, any- 
thing that touches it will be transformed into strange matter. Thus at first sight 
strange stars must have a bare quark surface ! This has to be a very sharply defined 
surface since it is limited by strong forces : at most it is a few fermis thick. More- 
over, due to the extremely high density of the surface (several times 10^'^ gm/cm^) 
the electron plasma frequency is so high that thermal radiation cannot be emitted. 
A bare strange star is a perfect silver sphere, not a black body. 

However, this may not be true. The quarks are confined by strong forces, but not 
the electrons. The electron sea thus extends slightly outside the quark matter, pro- 
ducing a charge disequilibrium at the surface which generates an enormous electric 
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field. Alcock, Farhi & Olinto estimate it to be around 10^'' — lO^^V/ cm. Positively 
charged ions can sit on top of the quark surface, suspended in the electric field, and 
form a crust of normal matter hiding the strange matter. The thickness of this crust 
depends on the strength of the electric field, i.e. on the electron density at the quark 
surface, which, in Farhi & JafFe's model of strange matter, in turn depends on the 
value of the parameters i?, and etc- There is however a maximum density that 
this crust can reach, independent of the model : as soon as neutrons drip out of the 
nuclei they will drop into the strange matter and be swallowed. Consequently, the 
crust maximum density must be below pdrip '■ a strange star can have a crust, but 
only an outer crust, no inner crust. This is the crucial difference between a strange 
star and a neutron star. 

It is not clear how a crust can form on a strange star. Matter falling radially 
on a bare quark star will have enough kinetic energy to reach the quark surface, 
but matter coming from an accretion disk could be stopped by the electric field. 
If a bare strange star is formed in a supernova explosion it is hard to believe that 
no material will fall back on it. However, once a very thin crust is formed, it will 
emit radiation probably at supra-Eddington rate and may blow out the crust. This 
certainly deserves more work and one cannot decide presently whether strange stars 
have a crust or not. 

We will assume in the next section that the strange star we study has a crust, 
of maximum thickness, i.e. with a density reaching pdrip- The mass of such a crust 
however is still very small, much less than lO""* Mq. 

4 Cooling of Strange Stars with a Crust 

We consider here the thermal evolution of a strange star with a crust and compare 
it to that of a neutron star. We assume it has been formed hot (supernova) and 
then cools freely (no reheating from accretion). Its thermal energy is lost by photon 
emission at the surface and by neutrino emission from the interior. For stars that 
are not too old the neutrino emission exceeds the photon emission by orders of 
magnitude. Fxcept when the star is hot (during the first few seconds of its life) the 
neutrinos produced have a mean free path much larger than the star radius and thus 
leave the star without any further interaction. Many different neutrino processes 
can occur in the core of a neutron star or in the strange matter, the most important 
ones being beta decay with inverse beta decay (URCA process) and variations of 
it. Table 1 shows the most common processes and their efficiencies : process a) 
occurs in all neutron stars, the occurence of processes b), c) and d) depends on the 
model. Process e) of course happens only in quark matter. These emissivities differ 
by several orders of magnitude and the quark processes are neither more nor less 
efficient than the nucleon ones, showing that the core neutrino emission cannot allow 
us to distinguish strange stars from neutron stars. Neutrino emission also occurs 
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Table 1: Some core neutrino emission processes and their emissivities (in 
erg / sec/ cm^). Ye = Ue/uB is the electron fraction and Tq the temperature in units 
of 10^ kelvins. 





Process Name 


Process 


Emissivity 


ci J 


IViOQineQ Ulvv^rV 


n + n' — 5- n' + p + e~ + z7^ 
n' -\- p -\- e~ — 5- n' + n + z/e 


^ 1 (-120 . 7^8 


b) 


K-condensate 
- URCA 


n + K~ — ^ n + e~ + z7^ 
n + e~ — ^ n + K~ + Ue 


^ 1024 . 7| 


c) 


TT - condensate 
- URCA 


n + 7r~ — 5-n + e~+z7^ 
n + e~ — 5- n + 7r~ + z/e 


^ 1026 . 7| 


d) 


Direct URCA 


n ^ p + e~ +17; 
p + e~ ^ n + Ue 


~ 1027 • r| 


e) 


Quark URCA 


J — 5- M + e~ + z7^ 

M + e~ — 5- J + Z/e 


~ 1027a,i;i/3r6 



in the crust but usually at a less efficient rate making that during the early phase 
of the cooling, the core is cooler than the crust. The surface is also much cooler 
than the crust due to the poor thermal conductivity of the material below it in the 
region where the electrons are not yet degenerate. Brown et. al.^ warn that the fast 
cooling of the core is not seen at the surface until the heat of the crust has had time 
to diffuse into the core and to the surface. From a simple random walk argument 
they showed that this diffusion time 

Td oc (crust thickness)2 (8) 

and is at least several decades for a neutron star. From this one immediately 
sees the difference between a neutron star and a strange star : since its crust is much 
thinner the diffusion time will be a few years for a strange star. Alcock & 01into2 
were the first to point this out. This was then confirmed by P. Pizzochero^ who 
recently showed how the heat diffusion in the crust can be neatly solved analytically, 
at the price of some simplifications in the involved physics. 

^.i Detailed Calculation of Strange Star Cooling 

To draw more quantitative conclusions we have studied in detail the heat prop- 
agation inside the star using a full stellar evolution code. The basic, general rela- 
tivistic equations are presented in the review by S. Tsuruta.'' The code used is an 
adaptation of the one developed by the author for neutrons stars, which includes 
the best crust physics (thermal conductivity, radiative opacity, neutrino emissivity. 
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specific lieat and equation of state) to date. For tlie quark; core, tlie foiiowing pliysics 
was used : 

• Tlie quark; E.O.S. is calculated following Farlii & JafFe"^ and Alcock, Farlii & 
Olinto'* by the method presented in sections 2 and 3. We chose twelve different 
sets of parameters (5, and etc) which span the whole range of values over 
which strange matter is stable. They are shown in Table 2 together with the 
properties of the corresponding stars. The crucial quantity obtained here is 
the electron concentration Ye because the neutrino emissivity is proportional 
to Yj-^^. All stars considered have a mass of 1.4 Mq. 

• The specific heat is simply given by relativistic Fermi liquid theory^'^ as (for 
each quark fiavor i) : 

Cv = 3vr^iV(0) ■T=^, pjii) (l + Aa.) • T (9) 

• The neutrino emissivity [see e) in Table 1] is from Iwamoto's work.^^ The 
reaction analog to e) with s quark instead of d is also included, but it is only 
a small correction since it is Cabibbo suppressed. 

• The thermal conductivity has been recently evaluated by Haensel & Jerzak,^^ 
we use their approximate formula : 

A = 3.4.10-(^)-'''(^)i ..,c™-„-A- ,10) 



^.^ Results 

The cooling curves of the twelve models are shown (as continuous lines) in Figure 
1. One sees that during the first months the surface temperature is the same for all 
models, this is because during this phase it is completely determined by the upper 
layers of the outer crust, which are the same for all models. After this phase the 
temperature drops, corresponding to the cold wave of the core reaching the surface. 
The time of this phase depends on the thickness of the crust, the thicker the crust the 
later the temperature drop (see Fq. 8). Figure 2 shows this explicitly with models 6 
and 12 : they have the same core cooling rate but different crust thicknesses as can 
be seen from Table 2. At the end of this period the crust and core are isothermal 
and we only have a temperature gradient at the surface. The temperature in this 
phase is entirely determined by the core physics and varies greatly from one model 
to another (see differences in Ye in Table 2). Later, photon cooling eventually takes 
over and the temperature differences disappear (this is not shown in the figure). 
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Table 2: Strange Stars (1.4 Mq) Models and Strange Matter Parameters. Rstar and 
Rcore are the star and core radii. During the coohng the envelope contracts shghtly, 
changing the star radius, but the core radius does not vary, pc is the central density. 
Ye = Ue/uB is the electron fraction, given here at the center of the star. 



# 


etc 

\C> CI / 


51/4 
Mev 


Mev 


Rstar 

km 


Rcore 

km 


P 




Ye X 10^ 


1 


0.3 


138 


100 


11.899 


11.461 


6.06 • 


10^4 


0.18 


2 


0.3 


138 


200 


11.444 


11.047 


7.18 • 


10^4 


8.2 


3 


0.3 


138 


300 


11.493 


11.091 


7.49 • 


10^4 


79. 


4 


0.3 


145 


100 


10.999 


10.640 


7.96 • 


10^4 


0.09 


5 


0.3 


145 


200 


10.553 


10.229 


9.63 • 


10^4 


4.5 


6 


0.3 


150 


100 


10.404 


10.092 


9.76 • 


10^4 


0.06 


7 


0.6 


128 


100 


13.327 


12.750 


4.16- 


10^4 


0.02 


8 


0.6 


128 


200 


12.761 


12.241 


4.96 • 


10^4 


2. 


9 


0.6 


128 


300 


12.707 


12.193 


5.40 • 


10^4 


35. 


10 


0.6 


135 


150 


12.013 


11.565 


5.93 • 


10^4 


0.08 


11 


0.6 


135 


200 


11.757 


11.332 


6.52 • 


10^4 


0.6 


12 


0.6 


140 


100 


11.621 


11.208 


6.55 • 


10^4 


0.003 




-101234 -101234 



Log t (years) Log t (years) 

Figure 1: Cooling of Strange Stars vs. Neutron Stars. The continuous curves 
correspond to the twelve models of strange stars from Table 2, the dashed curves 
to 1.4 M0neutron stars cooling with the neutrino processes a),b),c) and d) listed in 
Table 1. 

Figure 2: Flfect of Crust Thickness on the Cooling : comparison of models 6 and 
12. 
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The dashed hnes in Figure 1 are typical coohng curves of neutron stars with the 
various processes of Table 1. One sees that in the isothermal phase neutron stars can 
have surface temperatures which are comparable to strange stars and may not be 
distinguishable from them. However, all models of neutron star cooling (see®'^'^"^'^'* 
for more models) give a surface temperature above 2 millions kelvins for at least 
the first decade of the star's life, while all models of strange stars presented here 
result in a surface temperature below one million kelvins during the same period. 
This difference is a characteristic signature of strange stars, due to the thinness of 
their crusts and the fast cooling of their cores. The compact object recently formed 
by the core collapse in SN1987A is precisely of this critical age, raising the exciting 
possibility that detection of thermal radiation from its surface in the near future may 
unveil the presence of a strange star. 

One must also mention that the Vela pulsar has a temperature estimated^^ of 
just under one million kelvins for an age of about 10"* years, and thus it is well 
above the cooling curves of strange stars. Moreover A. Alpar^^ warned that strange 
star models have difficulties in explaining pulsar glitches, observed in particular on 
the Vela pulsar, providing a second argument for this object not to be a strange 
star. This does not disprove the strange matter hypothesis since the uncertainty 
about the possible mechanism which transforms of a neutron star into a strange 
star perhaps indicates that some of them do not metamorphosize. On the other 
hand, quark may be superfiuid,^'' and a mechanism similar to the one proposed by 
the author for kaon condensation®'^ may act here to keep strange stars warmer and 
could resolve the disprepancy with the Vela pusar. This aspect obviously requires 
more consideration. 



5 Conclusion 



The strange matter hypothesis of E. Witten claims that strange quark matter is 
the ground state of QCD at finite baryon number. If this were true it would imply 
that probably all neutron stars are actually strange stars. The possible masses and 
radii of strange stars are similar to those of neutron stars, as well as their rotational 
properties, making them impossible to distinguish from this point of view. 

We showed however that, due to the absence of an inner crust, strange stars with 
a crust have a very particular cooling history during the first decade(s) of their life 
and have then a much lower surface temperature than neutron stars. In reference 
to the SN1987A, the detection of thermal radiation from the surface of the compact 
object then formed, with a temperature below one million kelvins would be a definite 
signature of a strange star. 
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